clc, clear
data = importdata('../data/P1Element_64.txt');
[tmp, ~] = size(data);
n = sqrt(tmp);
x = data(:, 1); X1 = x(1 : n : end, 1);
y = data(:, 2); Y1 = y(1 : n, 1);
z = data(:, 3);
for i =  1 : n
   X(:, i) = X1(i, 1) * ones(n, 1);
   Y(i, :) = Y1(i, 1) * ones(1, n);
   Z(:, i) = z((i-1) * n + 1 : i * n, 1);
end
mesh(X, Y, Z)
xlabel('x');
ylabel('y');
zlabel('z');